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1. Introduction shown no interest in them because they cannot generate 


The focus on renewable energy sources has increased 
significantly in recent years due to environmental pollution, 
increasing energy demand, and declining fossil fuel resources. 
Various sources of renewable energy include biomass, solar, 
geothermal, hydroelectric, and wind. In the meantime, the 
wind has been proven to be a cheaper alternative to energy, 
and hence extensive research efforts have been made to 
improve wind power generation technology. The world has 
enormous potential for wind energy that can be used to 
generate electricity. According to research by Iran's New 
Energies Organization, Iran has areas with good potential for 
wind energy. Since wind speed is different in different parts 
of the country, to efficiently use wind energy, a suitable 
turbine must be selected with environmental conditions to 
have a proper economic justification [1]. Over the past 10 
years, electricity generation capacity has increased from 25 
GW to 200 GW [2]. The first type of man-made wind turbine 
is a type of drag-type vertical axis wind turbine (VAWT) used 
in the Sistan region [3]. Although VAWTs were the first 
turbines used by humans, modern-day researchers have 


electricity on a large scale. In recent decades, horizontal axis 
wind turbines (HAWT)s have been extensively researched, 
and a large part of the installed capacity of wind turbines is 
related to this type of turbine [4]. However, research on 
VAWT continued in parallel on a smaller scale. Solving the 
main problems related to VAWTs, such as blade lift force, 
blade failure due to oscillating force and oscillating wind, etc 
[5, 6]. are some common previous studies, but more research 
is needed. VAWTs are currently not economically attractive 
on a large scale [3]. However, VAWT is a way to generate 
energy in locations away from the main distribution line or 
where wind farms cannot be used due to environmental 
issues, so small-scale dispersed generation units are 
preferred. For this reason, mass production of VAWTs has 
only recently begun as a small-scale production unit [7]. 
Scientists have created various configurations of this type of 
turbine and have used different methods for their studies. In 
this paper, the details of these techniques and configurations 
are examined, along with the important findings of 
researchers about VAWTs, and a new method for 
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Nomenclature: 
A projected frontal area of turbine (m?) 
Ca blade drag coefficient 
Cop rotor drag coefficient Fp /pAV.2 
C blade lift coefficient 
F; tangential force 
Fp turbine overall drag force 


Fia average tangential force 

F* reduced frequency 

V speed (m/s) 

V; chordal velocity component 

Vou induced velocity in the upstream side 
u induction coefficient in the upstream 


Greek signs 

a AOA (rad) 

w angular velocity (rad/s) 

A= Ro /V,, tip speed ratio 

Subscript 

a induction 

d downstream condition 
Abbreviations 

HAWT horizontal axis wind turbine 

AOA angle of attack 

MST multiple streamtube 

BEM momentum models 

Ko-SST Kw-shear stress transport 

URANS unsteady Reynolds-averaged Navier-stokes 

DES detached eddy simulation 

SMM sliding Mesh Method 


C blade chord (N/m?) 

Cp turbine overall rotor drag coefficient 
Ch normal force coefficient 

Cp Power factor 

F, normal force (in the radial direction) 
H/C aspect ratio 

H height of turbine 

P Power (w) 

W Relative speed (m/s) 

Va normal velocity component 

Vaad induced velocity in the downstream side 
u' induction coefficient in the downstream 
p density (kg/m?) 

o solidity 

gen generator 

u upstream condition 

VAWT vertical axis wind turbine 

TSR tip speed ratio 

DMST double multiple streamtube 

PIV particle image velocimetry 

CFD computational fluid dynamics 

LES large eddy simulation 

RFM reference frame method 

CFL Courant_Friedrichs-Lewy number 


investigating this type of turbine will be proposed and 
compared with previous samples. 


2. Wind turbines configuration 

Differences in the shape and efficiency of wind turbines 
cause differences in the appropriate conditions for their use. 
In this section, we will review the types of turbines and 
compare them. 


2.1 HAWTs 

A group of turbines whose axis of rotation is horizontal. 
Early types include German and American windmills [6]. The 
first wind turbine in the United States was commercially used 
in 1931 by Charles Brush to generate electricity [5]. But 
modern VAWTs (Figure 1), which are used commercially for 
wind farms to generate electricity, usually have three blades 
that provide the force needed to rotate the rotor by lift force. 
The blades are positioned in the wind direction by a 
computer-controlled motor. At the top of the tower, there is a 
gearbox and a generator, and the blades are connected to this 
part. These turbines have been a major part of the 
contemporary history of wind energy technology and are a 
common commercial type. Although the efficiency of these 
turbines is high, it also has disadvantages that will be 
mentioned in the next section [8]. 


2.2 Drag-type VAWTs 

Sistan wind turbine, which is also the first type of wind 
turbine, falls into this category, but the most famous type of 
this category, which has a higher efficiency than other cases 
of drag-type, is the Savonius rotor, whose name is derived 
from the name of its Finnish designer. 


This type of turbine, which is a drag-type and vertical axis, 
was first tested in 1925. This rotor was very popular in the 
1960s and 1970s due to its simple construction and 
production, high starting torque, and ability to accept wind 
force from all directions. However, it has some disadvantages, 
such as high weight, low efficiency (15-19 percent), high 
torque changes on the axis of rotation, and wind from low 
latitudes which made it out of the spotlight over time made [9, 
10]. The mechanism of this type of turbine is such that the 
wind drag force in the concave part of the turbine is more than 
the convex part, and therefore it will cause the rotor to rotate 
[11]. After passing through the concave part, the fluid enters 
the other half of the rotor to apply a force against the direction 
of wind flow and in the direction of rotation of the rotor. 
Regardless of the drag force in the convex part, the turbine 
output capacity is calculated as Eq. 1 [5, 10, 12]. 


1 
P= Fau = 5 PAWair = Uturbine)” CpUturbine (1) 


And maximum power factor in the blade tip speed ratio (TSR) 
1/3: 


Cp =r = Cy (2) 


~ Spav3 27 
By comparing the maximum output power of the drag-type 
turbine with the Betz limit, the inherent weakness of the 
efficiency of the drag-type turbines is simply revealed. Even if 
we use the maximum drag coefficient, which belongs to the 
wind-facing hemisphere (Cp = 1.3), the maximum output 
power coefficient is 0.193, which is much lower than the Betz 
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limit. Another disadvantage of drag turbines is the efficient 
use of half of the space occupied by the rotor. Although wind 
technology began with drag-type mills, in practice, these 
inherent flaws have hindered the development of more 
advanced species [13, 14]. 


Savonius-Rotor Darrisus-Rotor H-Rotor HAWT 


*Meseuser™ 


ta, 


Figure 1. Different types of conventional commercial VAWTs 
[18] 


2.3 Lift-type VAWTs 

The turbine was invented in 1931 by Jean Marie 
Darrieus, a French engineer [6]. This turbine is of lift type, 
which means that, unlike the Savonius turbine, the lift force 
causes the turbine to rotate, and the drag force has a deterrent 
effect (Figure 1). The rotation speed of the Darrieus turbine is 
many times faster than the wind speed it blows. Therefore, 
the output torque of this turbine is low, and its rotation speed 
is high. Theoretically, this turbine has the same efficiency as 
the usual and commercial type of wind turbines, which in 
practice, considering the design limitations and resistance of 
the structure, it is difficult to achieve its efficiency. Using the 
lift force, the blade can move faster than the wind, thus 
producing more power [15]. The ratio of the blade tip to wind 
speed is a comparative indicator for wind turbines. This ratio 
is approximately 7 for turbine blades with maximum lift 
power, while 0.3 for blades with pressure drag [16, 17]. 
Figure 2 shows a comparison between turbine efficiencies. 
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Figure 2. Comparison of the power factor of different wind 
turbine 


When the airfoil moves in a circular path, it changes the 
relative speed of the wind relative to the airfoil, which 
ultimately produces a variable angle of attack (AOA). This 
wind movement on the airfoil produces lift force and, thus 
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torque. When the airfoil rotates around the main axis and 
passes through the opposite point, the AOA of the airfoil 
becomes negative, and thus the airfoil production force 
continues to produce torque by the direction of rotation [19, 
20]. In this rotor, the AOA with rotation is constantly 
changing. Therefore, in the front and back of the device, 
according to the wind direction, the maximum AOA and 
consequently, the most torque is produced. These pulsed 
torques make the design process complicated, and at some 
rotational speeds, the natural frequency of the blades equals 
this frequency, and this causes resonance and destruction of 
the device. Therefore, in design, these modes should be 
considered [19, 20]. Unlike HAWTs, where most of the weight 
of the blade structure is in the wheel bowl, in this rotor, most 
weight is located at the end of the arms [21]. Therefore, in this 
mechanism, large centrifugal forces are generated, and 
consequently, many stresses are applied to the structure, 
which adds to the complexity of the design. One of the 
initiatives to reduce this force is the use of rotating wings 
similar to the weight of an egg (also named troposkein) so 
that the blades themselves are attached directly to the base 
[22]. The rotating blades compress the column, and the 
column is secured from above with retaining wires. Although 
the Darrieus rotor is far more practical than the Savonius 
rotor, it still lags far behind horizontal axis turbines in terms 
of benefits and costs. Table 1 compares the mentioned rotors 
[23]. 


Table 1. Comparison of advantages and disadvantages of 
turbines 


Cases HAWT Darrieus Savonius 
es ore complicated simple simple 
Rectifier needs needless needless 
mechanism 
Ability to change 
the angle of the yes yes no 
blade 
Largeisral Yes Yes No 
construction 
Rotor weight medium low high 
bearing 
Noise high low low 
Generator and top of the rommá ound 
gearbox location tower 8 8 
Auto start yes weak yes 
Foundation complicated simple simple 
Ideal power factor high high low 
Maintenance costs high low low 
Research done much little little 
Structure set complicated simple simple 


3. Aerodynamic solution methods 

Although the flat-blade Darrieus turbine is one of the 
simplest types of wind turbines, its aerodynamic study is a bit 
complicated. Before comparing the existing designs for 
aerodynamic study, the equations used in these studies are 
introduced. The flow velocity is not constant upstream and 
downstream (Figure 3). The component of cord velocity and 
normal velocity is obtained as follows [24, 25]. 
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V; = Rw + Vacosé (3) 
Va = Vasin (4) 


In which, V, is the axial induction velocity that passes through 
the rotor, w rotor angular velocity and a is the rotor 
placement angle. According to Figure 3, the AOA can be 
calculated from the following equation [26]. 


-1% = sin 0 
a = tan * 7 = tan * Raj, (5) 
E Te +cos 6 
“Nix 


Also, the relative velocity on the blade is obtained as follows. 


Z 
WY, n_n (E Yæ) + cos 6] + sin?6 (6) 
Voo Va Væ Voo Vo” Vo 


The direction of lift and drag force, its radial and angular 
components, and the force coefficients in the angular and 
radial direction shown in Figure 3 are calculated as follows 
from the lift and drag coefficients obtained from the 
laboratory result [27]. 


C, = C, sina — Ca cos a (7) 


F, = 5 C,pCHW? (8) 


Given that Eq. 6 is written in terms of position, then the mean 
angular force can be defined as follows. 


2 
Fea == fy" Fe(8) a9 (9) 
And the output power is obtained as follows. 


P = F,,RNw (10) 


VER ~ ge 


É 


180° 
Figure 3. Velocities and forces on the Darrieus rotor [27] 


In the past, several mathematical models based on 
several theories have been used to predict the performance 
and design of the Darrieus turbine in various studies. 
According to research conducted by Islam [28], the most and 
best-validated models are divided into three categories: 

e Momentum models (BEM) 

e Vortex model 

e Cascade model 

Of course, none of these three models will meet all the key 
performance metrics. 


3.1 Momentum models (BEM) 
These models have been introduced in different ways, 
but the basis of these models is that by placing aerodynamic 
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forces in the direction of flow with momentum changes, the 
induced speed during flow will be obtained. The most 
important problem of these models is that they do not give a 
good answer for high TSR and high solidity because solving 
momentum equations alone is not enough [28]. In 1974, 
Templin [29] introduced the single streamtube (SST) model, 
which was the first and simplest. This model was first 
registered for the windmill but later expanded to the Darrieus 
turbine. In this theory, itis assumed that the induced velocity 
is constant along with the rotor and is obtained from the 
equality of the drag force along with the flow by a momentum 
change in this direction [30]. This theory pays attention to the 
effect of airfoil failure [31]. Thus, the effect of geometric 
variables such as solidity and aspect ratio as well as the effect 
of zero lift-drag coefficients are considered qualitative 
features in this theory, but the sheer effects of wind in this 
model are not considered [28]. According to Gluert's actuator 
disk theory [28], the uniform velocity passing through the 
rotor is equal to the average velocity upstream and 
downstream (Figure 4a). All these calculations are done by 
this model for a blade that equates the chord length of this 
blade to the sum of the blade lengths. The drag force in the 
flow direction is the result of momentum change in this 
direction, and by defining the rotor drag coefficient in Eq. 10, 
we have: 


Fp = M(Voo a W) (11) 
4(Voo—Va) 

Cpp = ig (12) 

Va _ 1 

Veo (=i) (13) 


The total torque for the turbine will be obtained from Eq. 
4-9 by placing va in Eq. 13 during an iterative process. This 


model can calculate the overall efficiency for the low load 
mode, but according to Wilson's results, it always achieves 
more output power than the laboratory results because the 
velocity on the downstream and upstream are assumed to be 
the same. This difference in results increases with increasing 
solidity and increasing TSR [28]. In 1947 Wilson [32] 
introduced the multiple streamtube (MST) model. In this 
model, the volume swept by the turbine is divided into several 
streamtubes side by side. The aerodynamic dependence of the 
parallel streamtubes is shown in Figure 4b Previous 
equations are written for each streamtube. In this model, the 
changes in induced velocity in front of the disk are considered 
in both vertical and horizontal directions. Thus, this model 
also includes wind shear effects [28, 33]. Wilson 
recommended the equation for the calculation of induced 
velocity as Eq. 14, taking into account only the lift force. 


ve = 1— (Ess «sing) (14) 


In 1981, Paraschivoiu [34] introduced the double-multiple 
streamtube (DMST) model. In this model, calculations are 
obtained for the two upstream and downstream sections 
according to Figure 4c The flow of air passing through the 
tube passes through two actuator disks in a row. 


Vau = Ue (15) 
Ve = (2u — 1)Vo0 (16) 
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Figure 4. (a) SST model, (b) MST model, (c) DMST model [28] 


Vaa = ul(2u — 1) (17) 


In which, Vau and Vaq are induced velocity in upstream and 
downstream side and, V, in the wake velocity. Also, u is the 
induction coefficient in the upstream, and v' is the induction 
coefficient in the downstream and is always u’ <u. The 
turbine's overall drag force in the direction of the stream is 
obtained from the following equation. 


o pe WwW? cos0 sin 
Fp E AeA |Įcos8| A Cr a) do (18) 
This method was able to solve one of Wilson's problems, 
which was the equality of the induction velocities of the front 
and rear blades, but problems with convergence (especially 
downstream and high TSR) appeared [28]. The results of 
these methods showed (Figure 5 and Figure 6) as the TSR 
increases, the turbine power coefficient increases to reach an 
optimal TSR. After this, the power factor will decrease with 
increasing TSR. Paraschivoiu was also able to study the effect 
of solidity as shown in Figure 7. His results showed that with 
increasing solidity, the power coefficient in TSR reaches its 
optimal value. 
0.5 
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04 ==- Multiple Streamtube 
Experimental data 

a © Vo=7 m/s 
003) O Vo= 9 ms 
5 © Vo= 11 m/s 
o 
& ii 
2 02 Solidity = 0.27 
5 
5 
è 


Figure 5. Comparison of SST and MST models with laboratory 
information [35] 


February 2024] Volume 03 | Issue 01| Pages 23-33 


(c) 


In 2013, Batista et al. [35]. made a comparison between 
the models as shown in Figure 5 and Figure 6. His results 
showed that the DMST model has the closest adaptation of the 
three models to laboratory information but also, has the most 
difficult convergence conditions. All three models at higher 
solidity and higher TSR predict higher results from laboratory 
data. With the help of these methods, it was possible to obtain 
the main aerodynamic properties with low calculations, but 
these methods, which were not able to model the shape of the 
flow, are unable to investigate the causes of phenomena such 
as power factor reduction for TSR greater than 6. Another 
problem with these methods is the lack of proper airfoil 
information at attack angles above the stall state. 


0.6 
© SANDIA field test Prediction data (Re = 0.3x105) 
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Figure 6. Comparison of multiple and DMST models with 
laboratory information [35] 


3.2 Vortex model 

The vortex model is based on potential flow. This model 
calculates the velocity field around the turbine through the 
influence of velocity in the wake of the models. Turbine blades 
are expressed as concentrated vortexes. The strength of these 
vortexes is obtained by collecting airfoil coefficient data and 
calculating the relative flow and AOA. In VAWTs, the blade 
element is replaced by a substitution vortex filament or a lift- 
line [28]. This method is well compatible with laboratory 
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information and was able to solve some of the problems of the 
previous method. Dixon's results [37] could explain that as 
the blade speed increases, less energy is wasted from the 
wind, but after exceeding 6 times the wind speed, the blade 
collides with the non-renewed wind flow, which reduces the 
power factor (Figure 8). However, this method cannot 
visualize the flow well and requires a lot of calculations [28, 
38]. 
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Figure 7. Turbine power coefficient curve in terms of TSR for 
o=0.8 [36] 
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Figure 8. Displays the rise motion of the vortex method [37] 


4. Flow field visualization by particle image 
velocimetry (PIV) 

Aerodynamic models, no matter how well they predict 
the power and performance coefficients of the turbine, are 
still unable to visualize the flow field. Visualization of the flow 
field around the wind turbine facilitates the understanding 
and analysis of the aerodynamic behavior of the wind turbine, 
which is necessary to improve the efficiency and proper 
design of the device [39]. In 1999, Fujisawa et al. [40] 
visualized the flow field in an area of the Darrieus rotor 
during dynamic stalls by injecting dye. In this method, the dye 
was injected into the fluid stream, and the flow model was 
observed through the dye path. They used the PIV method and 
the regional imaging technique to measure the average 
velocity distribution around the blade. In PIV, tracer particles 
with a specific gravity close to the fluid flow are injected into 
the fluid without affecting its velocity. These particles move at 
the same local velocity of the fluid. Images of moving particles 
are taken at any given moment with the help of a digital 
camera, and then these images are analyzed to obtain the 
velocity of the fluid at that moment. They concluded that the 
dynamic stall phenomenon occurs due to the flow of two pairs 
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of vortices from the blades as they rotate from the rotor. At 
low rotor speeds (TSR<3), the AOA is variable and dynamic 
stall occurs [39]. In 2014, Danao et al. [41]. compared the 
turbulence models Ke and Kw-SST with laboratory results 
(Figure 9). Their results showed that model Kw-SST has a 
greater ability to model the flow around the rotor. Model Ke 
was initially introduced for high-spin currents, but over time, 
corrections were made to model the flow along the wall. 
Therefore, as previously predicted according to their results, 
this model is better consistent with the laboratory results 
around turbomachines and because it is two equations does 
not increase the calculations. 


O CDF:kw SST 


CDF:k-e standard PIV 


70° 60° 


140° 


Figure 9. Flow visualization using Ke, Kw-SST and PIV [41] 


In addition, they examined the effect of wind blowing on an 
unsteady basis. Their results showed that if the wind enters 
the turbine in an oscillating manner, the optimal of the blades 
differs and occurs in less than the steady-state (Figure10). 


0.5 —- -=-= 


steady CP 


CP 


Figure 10. Power factor curve for both steady and oscillating 
wind modes [41] 


5. Computational fluid dynamics (CFD) methods 

CFD is now becoming a powerful tool in fluid mechanics 
that analyzes and solves fluid flow problems using numerical 
methods and algorithms with the help of electronic 
computers [42]. Using CFDs can save time and costly testing 
and can also be used to improve VAWT analysis. A well-tuned 
CFD model can simulate actual flow conditions that can 
produce results consistent with laboratory outputs. 
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Therefore, design optimization can be done before entering 
the construction stage [43]. CFD analysis has been widely 
used to determine power factor in different types of wind 
turbine configurations. The CFD technique for visualizing the 
flow field of a VAWT can be seen in Allet's early work [44]. 
Debnath et al. [45] used the CFD to test the Savonius-Darrieus 
rotor, which is more capable of starting than the Darrieus 
rotor. In 2010, Wang et al. [46] used two unsteady Reynolds- 
averaged Navier-stokes (URANS) models with standard Kw 
and Kw-SST to simulate the phenomenon of the dynamic stall 
at low Reynolds (Rex105). After comparing the results of this 
model with the experimental model, they found that model 
Ko-SST is better than standard Kw. Model Kw-SST shows the 
main features of the dynamic stall phenomenon, such as the 
aerodynamic residual load and the flow structure of the 
vortexes of the front of the airfoil. Large eddy simulation 
(LES) models are more advanced and more powerful than 
URANS models. They require more computation time but 
their results are more accurate. Among the two models, 
detached eddy simulation (DES) and LES, the DES model is a 
better model due to the greater agreement of the results with 
the experimental results. The DES model is a combination of 
the LES model and the URANS model. This model not only has 
a lower computational cost than the LES model but also 
models the wall area more accurately [47-49]. To test 
turbines by the CFD method, in addition to the equations 
governing the fluid flow field, the movement of the turbine 
blades must also be considered. In general, turbine motion is 
modeled in three ways: stationary rotor, constant speed 
rotor, and accelerated rotor [50]. 


5.1 Stationary rotor 

In this method, it is assumed that the rotor is fixed and 
the resulting current and torque are obtained for several 
angles of the fluid field. Menet et al. [51] compared the flow 
field for several different angles of the two Savonius turbines 
but did not obtain good values for the power factor compared 
to the laboratory data. Of course, this method is not very 
useful for lift-type turbines such as Darrieus due to the high 
importance of the relative speed on the rotor. In 1999, Allet et 
al. [44] investigated the behavior of an airfoil at higher angles 
of the stall for which it did not have adequate laboratory 
information at rest and combined aerodynamic methods to 
obtain the power factor for Darrieus turbines. Of course, their 
data always showed a lower power factor than laboratory 
data, due to the inability of this method to consider dynamic 
stall resulting from its static study [50]. 


5.2 Constant speed rotor 

In general, three methods of reference frame method 
(RFM) and sliding mesh method (SMM), and dynamic mesh 
are often used to study turbomachines that have moving parts 
[52]. In the RFM method, which is also the simplest method, 
the relative coordinate system is adapted to the rotor, and 
instead of the mesh deformation, it is assumed that the 
boundary conditions are rotating around the axis and the 
momentum equations for the rotating device are correct. Of 
course, this method is not very useful for analyzing 
turbomachines because it solves the problem steadily and 
does not give a good analysis of the problem, which is 
generally transient [8, 53, 54]. Korobenko et al. [55] used a 
complete method known as SMM. In this method, unlike RFM, 
the mesh shape will change. Amet et al. Divided the solution 
area into a fixed part and a rotating part in which the rotor 
was located. For each subdomain, the governing equations 
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are solved and related by Eq. 19 and 20 through our boundary 
which is known as the interface boundary. 


Uy — Us = 0 (19) 


In which, I is the unit tensor, S and M are the fixed amplitude 
and the moving amplitude, respectively, and n are the unit 
vectors perpendicular to the outside. In this solution, the 
rotational amplitude equations are given a stationary angular 
velocity and the forces and torque on the turbine are 
calculated at each time step. In 2009, Amet et al. [56] used the 
same method to investigate the effect of a vortex on the 
operation of a Darrieus turbine. By defining the reduced 
frequency as Eq. 21, they introduced the ratio of the time of 
fluid movement on the blade to the time scale of the intensity 
of changes in the AOA. 


(21) 
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Figure 11. a curve in terms of 9 for different A and F* in terms 
of A for different solidities [56] 


For high F* (greater than 0.05), if the AOA is greater than 
the static state, the device is prone to strong dynamic failure 
with vortex separation. Figure 11 shows two curves, first a in 
6 for different A and the other F* in A for different solidities. 
The first curve shows that 2<A<6, results in the 9.5<AOA<30 
degrees. These angles are higher than the 12 for the static 
stall. For A>5, the AOA does not exceed 12 in the whole half. 
The second curve shows that if the two conditions A>5 and 
C/R>0.1 are met, a strong dynamic failure is expected. Amet 
compared the separated vortexes for the two AOA 2 and 7. 
Figure 12 shows that the pair of vortexes that are clockwise 
and counter-clockwise detached from the blade when it is 
upstream, and in the downstream part of the rotor may 
collide with the blade, causing a drop in the blade lift 
coefficient. Of course, whether this reduction in lift factor 
reduces the rotor power depends on whether the torque in 
that area is positive or negative. Figure 13 shown for A=2, 
shows an oscillation drop-in area 8 due to the movement ofa 
vortex along the stream, which continues in area 15. Figure 
14 shows that if the TSR increases, fewer vortexes are formed 
with lower power and have little effect on the lower blade. 
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Figure 15 shows the lift factor in a complete rotor rotation in 
A=2. According to this diagram, with the increase of a, the 
oscillation caused by the collision of vortexes with the blade 
is no longer observed. Amet also showed that by increasing 
the TSR, the flow is not able to repress the pressure, and the 
turbine efficiency decreases [56]. 
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Figure 13. Lift coefficient in a complete revolution, in A=2 
[56] 
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Figure 14. The path of separated vortexes in A=7 [56] 
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Figure 15. Lift coefficient in a complete revolution, in A=7 
[56] 


6. One degree of freedom rotor 

According to studies, all the research that has been done 
so far has been done with the constant rotation speed method, 
and only Alessandro in 2010 [57] and Jaohindy in 2013 [58] 
used the accelerated rotation speed for the Darrieus rotor. In 
this method, the torque applied to the turbine causes 
proportional acceleration and motion in the turbine, which 
means that the rotor motion will be obtained by coupling the 
Lagrangian equations of CFD and the Euler equations of the 
rigid body. Figure 16 shows the dynamic solution flowchart, 
which consists of three main processes; Solve fluid flow, rotor 
processing, and mesh deformation processing. In other 
words, according to this method, the motion of a solid body is 
coupled by solving the flow field. According to the flowchart, 
the initial angular velocity is first considered zero. Then the 
equations governing the fluid field are solved, and the forces 
acting on the rotor are calculated. In the next step, the angular 
acceleration will be calculated by the code written according 
to Eq. 22. 


Joz8 + Moz (22) 
Moz = Mp +M, + M; (23) 
M; = Lk=0 a,w* + Mgen (24) 


In which, J,, is the moment of inertia of the rotor around the 
z-axis, M, is the torque due to pressure, M; is the torque due 
to viscous forces and, Mr, is the torque generated by the 
generator load, which also includes mechanical friction. It 
should be noted that M, and M, are obtained by calculations 
from the fluid flow field calculated in the previous step, and 
M, is obtained by Eq. 24, whose coefficients are taken into 
account by the laboratory information. Using the acceleration 
obtained in this step, the new location of the rotor is 
determined in the next time step. In the next step, the 
appropriate mesh is deformed with the new location of the 
rotor by the SMM method [58]. Since the time step is a 
function of the Courant_Friedrichs-Lewy number (CFL) 
number and the rotation speed of the rotor, it is 
recommended to change the time step according to the rotor 
for the appropriate computational volume. Fluent software 
can calculate the time step in each time step by writing code 
[53, 59, 60]. 
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As mentioned, the main weakness of the Darrieus rotor is its 
initial start. DeCost et al. [61] by examining the start of the 
Darrieus rotor by the PID method found that the drag forces 
had a greater effect on the start of the rotor than the lift force, 
which may even rotate in the opposite direction, but after 
increasing speed and the effects of lift force, the rotor 
accelerates in the right direction to reach its working speed. 
Flow field studies can be useful to obtain optimal starting 
conditions. Conventional methods that keep the rotor speed 
constant can certainly not work for this mode, and there is a 
need for a method that can study the rotor speed based on 
actual conditions. In addition to the starting period of the 
turbine in the rotor operating mode, the speed is not constant 
even during one rotation. Based on the results of Jaohindy et 
al. [58] it can be said that after passing the starting period 
during a rotation, due to the change in torque at different 
angles, different accelerations are obtained from Eq. 22, and 
as a result, we have different angular velocities during a 
rotation. According to Eq. 6, angular velocity strongly affects 
the relative velocity and lift force, which results in torque with 
an oscillating amplitude greater than the real state because 
the rotor speed is obtained from the flow conditions. 
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Figure 16. The flowchart is proposed for the dynamic 
solution method [58] 


7. Conclusion 

In this article, wind turbines were first compared. 
Darrieus turbine, along with many advantages, has a problem 
of starting and efficiency less than HAWTs, which can be 
widely used if these problems are solved. The various 
methods used to evaluate the aerodynamic characteristics of 
VAWTSs were reviewed. Meanwhile, momentum models are 
cheap but inaccurate and are considered suitable for iterative 
optimization algorithms. The Vertex model has acceptable 
accuracy but has a high computational cost, which has 
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prevented the use of these models. Of course, with the 
advances made, these models also have some ability to 
imagine the flow. Many CFD studies have been performed on 
Darrieus turbines that have been able to visualize the flow 
field well to understand the aerodynamic phenomena around 
these turbines. Understanding the flow is necessary to 
optimize the turbine, but most researches have assumed that 
the rotor speed is constant. In general rotor researches, this 
simplification does not cause a significant error, but these 
models are incapable if special conditions such as the start 
time of the rotor should be examined. Finally, it is suggested 
that according to an iterative algorithm, the rotor velocity be 
obtained as the actual torque reaction. 
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